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^ ■ The theory of the on-shell Sudakov form factor to all order of 
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^ I The key to understanding and using perturbative QCD is the idea of factorization. 
p i' Factorization is the property that some cross-section or amplitude is a product 
p • of two (or more) factors and that each factor depends only on physics happening 
1^ . on one momentum (or distance) scale. The process is supposed to involve some 
large momentum transfer, on a scale Q, and corrections to the factorized form are 
suppressed by a power of Q. (In general the product is in the sense of a matrix 
^ ' product or of a convolution.) 

The standard factorization theorems are typified by the one for the moments 
of the deep inelastic structure functions: 



1 INTRODUCTION 



FniQ) = Cn («s(Q)) * exp 



, 7„ (as(/u)) d/i//i 
Qo 



*M„(Qo). (1) 



Here F„ is the nth moment of one of the structure functions, C„ is a Wilson coeffi- 
cient, Mn is a hadronic matrix element of an operator of spin n and twist 2, 7„ is 
an anomalous dimension, and Qo is a fixed scale. The symbol denotes a matrix 
product (to allow the possibility of contributions from more than one operator). 
The renormalization group has been used to absorb all logarithms of large mass ra- 
tios into the integral over the anomalous dimension. These theorems are described 
elsewhere in this volume. 
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In this article, I will treat the Sudakov form factor, which provides the sim- 
plest example of factorization theorems of a more complicated kind. The difference 
between this case and deep inelastic scattering results from a difference in the regions 
of loop-momentum space that give the leading-twist contributions to the process. 
In the case of simple factorization theorems like (1), these regions involve lines with 
momenta that are either collinear to the detected particles or are far off-shell. In 
individual graphs there are leadin g tw ist contributions from regions with soft gluons; 
but after an intricate cancellatiorM^^, the effects of soft gluons cancel. However, in 
the Sudakov form factor the effects of the soft gluons do not cancel. Even so, a more 
general factorization theorem holds for this case. 

This form factor is the elastic form factor of an elementary particle in an 
abelian gauge theory at large momentum transfer Q. SudakoJ^ treated the off-shell 
form factor in the leading logarithmic approximation. I will treat the on-shell case 
and derive the full factorization formula, which is valid to all orders of perturbation 
theory and includes all nonleading logarithms. 

This relatively simple case is a prototype for such processes as the Drell-Yan 
cross section when the transverse momentum is much less than the invari ant mass 
of the Drell-Yan pair. A treatment to all order s of logarithms is given in 1^1^, and 
applications to phenomenology can be found in 1^1^. Work at the leading logarithm 
level can be found inl^^ and references therein. In^^ Sen showed how to treat the 
on-shell Sudakov form factor in a non-abelian theory to all orders of logarithms. 

The first step in proving any factorization theorem is to understand the re- 
gions of the space of loop momenta that give the "leading-twist" contributions, i.e. 
contributions not suppressed by a power of Q. After appropriate approximations it is 
possible to use Ward identities to convert the leading-twist contributions into a form 
that corresponds to the factorization theorem. A complication here is to eliminate 
double counting. Finally a differential equation for the evolution of the form factor 
is derived. It is only after this step that it is possible to perform systematic pertur- 
bative calculations, without having the validity of a finite-order calculation being 
brought into doubt by the possibility of large logarithmic corrections in higher order. 
The solution to the equation is in terms of quantities with perturbation expansions 
that have no large logarithmic terms in their coefficients. 

One topic I will emphasize is the extent to which the fact that we are dealing 
with a renormalizable gauge theory of physics comes into the form of the factor- 
ization. To do this I will start by examining what sort of result holds in a super- 
renormalizable theory without gauge fields, specifically 0^ theory in four space-time 
dimensions. The {(t)^)^ theory is of course completely unphysical. However, it is a 
simple model which exhibits the features common to any superrenormalizable the- 
ory without gauge fields but which has no irrelevant complications. We will see that 
a very simple factorization holds true. 
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Next, I will step up the space-time dimension to d = 6. This will render the 
model merely renormalizable instead of superrenormalizable. The short-distance 
part of the factorization will then become non-trivial, but it will still be of the same 
form as for deep inelastic scattering, (1). 

Finally, I will return to four dimensions, but now with a gauge theory. To 
provide a simple demonstration of how and why all the logarithms are under control, 
I must avoid complications that are irrelevant for this purpose. So I will take the 
theory to be abelian, with a massive gluon, and treat the annihilation of a qq pair 
into a virtual photon. The complications thereby avoided include: a non-abelian 
gauge group, the infra-red divergences caused by a massless gluon (which would 
mean we would have to discuss a kinematically more complicated process), and 
color confinement (which would force us to treat, say, a form factor of a composite 
particle). Although these complications are important for real strong interactions, 
they are inessential if we are trying to understand "Sudakov" effects by themselves. 



2 REDUCED GRAPHS 



Consider a form factor 

F={0\jmPAPB)- (2) 

Here j is a composite field, for example the electromagnetic current of quarks in 
QCD, and \paPb) represents an incoming two-particle state with energy Q: 

Q' = {pa+Pb)\ (3) 

which we assume to be very large. 

First we must find the regions of momentum space that are important in 
Fevnma n g raphs for this amplitude. We use the method given by Libby and 
StermarM-iJ. Suppose we scale all momenta by a factor Q: 

m = mQ. 

The reduced mass m goes to zero as Q goes to infinity, so that we are effectively going 
to a massless theory. If all scaled momenta in a graph are off-shell by order unity, 
then we get a contribution of order unity (given that we are in a renormalizable 
theory, so that the coupling is dimensionless) . Then simple perturbation theory is 
applicable provided only that the effective coupling as{Q) is small. Other leading 
contributions can come from regions where some of the scaled momenta become 
on-shell in the massless theory (m = 0). In this case we obtain a contribution only 
when the contour of integration is trapped at the on-shell point, for otherwise we may 
deform the contours into the off-shell region. Such points were called pinch-singular 
points inim. 
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Figure 1: Typical reduced graph for annihilation form factor. Jet lines are solid, 
and soft lines are dashed. 

The analysis of Coleman and NortorP^ can be used to locate the pinch 
singular points. Each pinch-singular point can be represented by a reduced graph. 
The lines whose scaled momenta are off-shell by order unity are all contracted to 
points; they form the vertices of the reduced graph for a pinch singular point. The 
lines that have on-shell scaled momenta form the lines of the reduced graph. In 
order that the contours of integration be pinched, the reduced graph must represent 
a classical scattering process. In the case of an annihilation form factor, the reduced 
graphs have the form exemplified by Fig. 1. The on-shell lines either have non-zero 
fractions of the scaled momentum of one or other of the incoming lines or they 
have zero scaled momentum. These lines are represented by solid and dashed lines 
(respectively), and are called jet and soft lines. An arbitrary number of jet lines 
parallel to interact and enter the reduced vertex where the annihilation occurs. 
A similar situation occurs for ps- An arbitrary number of soft lines join the two jet 
subgraphs. 

This can all be said without knowing the field theory. We next need to know 
which of the pinch singular points give important contributions as Q — > oo. For 
this purpose we consider only leading twist contributions, i.e., those that are not 
suppressed by a power of Q. Which regions give leading- twist contributions will 
depend on the theory within which we work, especially on its renormalizability or 
superrenormalizability and on the presence or absence of gauge particles. 
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(a) 



(b) 




(c) 



Figure 2: One loop graphs for form factor in 0^ theory. 
3 SUPERRENORMALIZABLE SCALAR THEORY: (cj)^)^ 



The Lagrangian of 0^ theory is 



(5) 



In d = 4 space-time dimensions the couphng, g, has positive mass dimension. This 
signals that no infinite coupling or wave function renormalization is needed, i.e., that 
the theory is supcrrcnormalizable. We define the form factor in eq. (2) by choosing 
the composite field j to be ^0^- 

First consider the one- loop graphs, which arc listed in Fig. 2. The graphs with 
self-energy corrections. Fig. 2(a) and (b), have reduced graphs equal to themselves. 
If these and higher-order self-energy graphs were all that we have, then the form 
factor would be equal to 

z{m'^,g), (6) 
where z is the residue of the renormalized propagator: 



Sf{p ,m ,g) 



p2 _ 



as p 



m 



(7) 



In fact these graphs are all that we have, for vertex graphs like Fig. 2(c) 
vanish as Q — > oo, by a power Q, as we will now show. Fig. 2(c) has the value 

= ^/ d^^ .„.. (8) 



(27r) 



-ig" 



(27r)4g2 



[w? — k"^) \w? — [pa + kY] [to? — {pb — kY] 

1 



d^k- 



k^) 



{pa -\-kY m? - {pb - ky 



(9) 
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The possible reduced graphs for Fig. 2(c) are hsted in Fig. 3: 

(a) Fig. 3(a) corresponds to the region that j/c'^l ~ 1, i.e. where all internal lines 
are far off-shell. Manifestly the resulting contribution is 0{1/Q^). This is 
essentially a result of dimensional analysis coupled with the positive dimension 
of g. 

(b) Fig. 3(b) corresponds to the region where is coUinear to p^. This graph is 
also O (1/Q^). To see this, we use light cone coordinates 

Pb-Pa-^^, Pb-P. = ^2- (10) 

Then in the region symbolized by Fig. 3(b), we have 

k+ = 0{Q), k-=0{X^Q), kT = 0{XQ), (11) 

where A is small. It is now easy to check that the contribution to the form 
factor is O {1/Q^). The point is that one quark line is far off-shell and that 
there are no compensating numerator factors. 

(c) Fig. 3(c) is just Fig. 3(b) with A ^ B. 

(d) For Fig. 3(d), which corresponds to the region where all components of k^ are 
much less than Q, we let 

P = O (AQ) , 

for all components. Again we get a contribution of order 1/Q^. 



It is fairly easy to show that this analysis holds true to all orders, by the 
methods of The reduced graphs for the leading twist contributions all have the 
form of self energy graphs attached to the lowest order vertex, so that the form 
factor is z plus higher-twist contributions. (See eq. (7) for the definition of z.) 



4 RENORMALIZABLE SCALAR THEORY: (0^)6 



4.1 One-loop 

The sole significant difference in going to six space-time dimensions is caused by the 
coupling's becoming dimensionless and the consequent need for coupling and wave 
function renormalization. We write the Lagrangian in the form: 



+ Uz{dct>f - \bm^^ - \^'^g0- + H. 
I Z b 



(12) 
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(d) 

Figure 3: Reduced graphs for Fig. 2(c). 

Here m and g are the renormahzed mass and couphng, and the last four terms 
are the renormahzation counterterms. We will use dimensional regularization (i.e., 
space-time dimension d = Q — 2e) to cut off the ultra-violet divergences. To 
keep the coupling g dimensionless we introduce the unit of mas /i. The lin- 
ear term hip is adjusted to cancel tadpole graphs; for the other terms we will use 
MS renormalizatiorP^. 

The structure of the reduced graphs is the same, as always. What changes is 
the size of the contributions. Consider the one- loop vertex graph Fig. 2(c), whose 
reduced graphs are in Fig. 3. By following the same method as we used in Sec. 3, 
we find that the contributions of the reduced graphs (at d = Q) are 

Fig. 3(a) Q°, 
Fig. 3(b) or (c) \^Q\ (13) 
Fig. 3(d) A^go. 

Clearly we get a leading contributions solely from the region where all internal lines 
are far off-shell. The existence of this contribution is tied to the dimensionlessness 
of the coupling. 

The contribution of Fig. 2(c) is therefore given by neglecting all masses, with 
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Figure 4: General graph for form factor. The cross-hatched bubble is the sum of all 
IPI graphs for the form factor. 



errors of order Thus 




+ counterterm. (14) 




(15) 



After using MS renormalization to cancel the ultraviolet divergence in the 
integral in (14), we find 



(At one-loop order, the MS scheme is defined by requiring counterterms to be a 
coefficient times 1/e — 7 + ln(47r), where 7 is Euler's constant. In the MS scheme 
we would omit the 7 and the ln(47r).) 

4.2 Higher orders 

The generalization to all orders of the one-loop results is obtained by observing that 
the graphs for the form factor are a product of two propagators and a one-particle- 




ln(-QV/i2) + 3 +0(l/g'). 



(16) 
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Figure 5: Leading reduced graphs for form factor in {4>^)q theory. 



irreducible (IPI) vertex (Fig. 4). A leading contribution is only obtained from the 
IPI vertex when all its internal lines are off-she ll by order Q^. It can be shown fairly 
easily that other regions are power suppressed^. Thus the only leading reduced 
graphs have the form of Fig. 5. This result is true in any renormalizable non- gauge 
theory. 

Therefore the form factor has the form 



Here F is the IPI vertex with the masses set to zero. This is the simplest example of 
a factorization theorem. The z factor comes from the single-particle propagator; it 
depends on phenomena on the scale of the quark mass m and is independent of the 
energy Q. On the other hand the vertex factor F depends only on the total energy 
and not on the mass. 

Since the theory needs renormalization, there is important dependence on 
the unit of mass fi. A perturbative calculation of 2; or F has large logarithms of 
m/ fi or of Q/ fi, so a simultaneous direct calculation of these quantities to low order 
cannot be reliable if Q/m is large enough. The renormalizatio n grou p comes to our 
aid since both z and F satisfy renormalization group equationJ^^^I^ 



F = z{g, m//i) F(g//i, g) + O ((logarithms of Q)/Q^) . 



(17) 



d 



27 {9 {^^)) 



(18) 




7<^2 {g (/i)) + 27 {g (/i)) 



7 id (/^))- 



(19) 
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Here 'j{g) and 7^2 (gf) are the anomalous dimensions of the operators (f) and 
spectively, and the renormahzation-group operator is 



re- 



d 



(20) 



dm? 



Note that 7 is half its value as defined by many authors. 

Evidently we can solve Eqs. (18) and (19) and write 



r /""^ du 

F(Q, m, g, /x) = z{g{m), 1) exp 2 / —7 {g{^')) 



ew\ ^%li9i^^')) r(i,^(g)). 



(21) 



Here g{fi) is the running (or effective) coupling at scale /i. Evidently each factor 
may be reliably calculated without large logarithms in higher-order corrections. 

4.3 "Optimization" of perturbation calculations 

In (21) the endpoints of the integrals over /i' are fj,' = m and /x' = Q. This is not 
necessary; all that is needed is that the endpoints be of order m and Q. This is 
important in "optimizing" perturbative calculations. We can write 



Here ci and C2 are arbitrary constants to be chosen at will. In a calculation to all 
orders of perturbation theory the result for F is independent of our choice of ci and 
C2. But in a finite order calculation the result has dependence on ci and C2 of the 
order of the first uncalculated term. We should choose Ci and C2 not too far from 
unity to keep higher order corrections small. The change in F given by varying ci 
and C2 by a factor of 2 can be regarded as an estimate of the error in the calculation 
induced by uncalculated higher order terms. 

There has been much discussioE^^ of appropriate ways to choose ci and C2. 




(22) 
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5 GAUGE THEORIES 



We now consider a form factor in the massive abelian gauge theory whose Lagrangian 
is 

£ = - -Fl, + \m^A^ - la . + g (z^ + ^^^g^ - M) q 

4 ^ (23) 

+ UV counterterms. 

Here the notation is standard. The renormahzed masses and couphng are m, M 
and g. We regulate ultra-violet divergences by continuing to space-time dimension 
(i = 4 — 2e, and we will call the gluon field and q the quark field. We will treat 
the electromagnetic form factor of the quark^ 

F = mqi^<imU{PA),q{PB)) , (24) 



when the center-of-mass energy Q [ = \l {pA+ PbY) gets large compared to all 
masses. 

SudakoJSl was the first to discuss such a form factor to all orders of pertur- 
bation theory. His result was for the sum of all the leading logarithms, but with 
the quarks off-shell. He found that F ~ exp —{g'^/STr'^) \n^{Q'^) . The on-shell case, 
with a massive photon, was first treated (still in leading logarithm approximation) 
by Jacki\\II^, with the result that F ~ exp l-i^g"^ /IGir'^) \n^{Q^) . 

From Sudakov's work (1956) until 1980, there was no progress in g oing sys- 
tematically beyond leading logarithms, despite many attempts. Muellep^ and 
CollinJSU' then gave an all-orders and all- logarithms treatment. The treatment be- 
low is an improved version of 1^ . A first version of the present treatment appearecP21 
as notes on lectures given in 1984. 

Notice that in these gauge theory form factors there are two logarithms of Q 
per loop rather than the one logarithm per loop that we have in {(f)^)^ theory. This 
is a symptom of the new physics present in a gauge theory. The effects that we will 
investigate reappear in many processes in QCD. 



5.1 One loop 

Self-energy graphs contribute just as they do in (p^ theory, and give an overall factor 
Z2{g,m, M, fi), the residue of the pole of the quark propagator. So the only non- 

^To be precise, note that the operator [q'y'^q] in (24) is the renormalized operator qj'^q + 
UV counterterms. 
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Figure 6: One loop vertex in QED. 



trivial one-loop graph is the vertex graph, Fig. 6. Its value is^ 




167r4 



(2vr/i) 



2e 




[-{PB - kf + M2 - ie] (m2 - k"^ - ie) [-{pA + kf + - le] 



(25) 



where the numerator of the gluon propagator is 



Nx,{k) = gxu - kxk,{l - 0/ {k^ - im^ + ^e) . 



(26) 



The possible reduced graphs are exactly the same as in 0^ theory and are listed in 



Since our theory has a dimensionless coupling, the ultraviolet region, of large 
fc, contributes to the leading power of Q, just as in 0^ theory in six dimensions. 
However, unlike 0^ theory in either dimension, the other three regions also give 
leading contributions. This happens because the numerator factor in (25) is 0{Q^) 
in all three regions; we know from our analysis of cj)^ theory that the graph would 
otherwise be of order We must first understand exactly how this factor of 

arises, since a systematic treatment of such enhancements by numerators is the key 
to a complete treatment of the form factor to all orders of perturbation theory. The 
mere fact that the reduced graphs Fig. 3(b), (c) and (d) are leading twist means 
that a simple factorization like (17) or (21) cannot hold. 

^Our 7-matrices are those of Bjorken and DrelPS, except that our wave functions satisfy the 
normaUzation conditions uu ~ 2m, vv — —2m. 



Fig. 3. 
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5.2 Method of Grammer and Yennie 



Consider the region corresponding to the reduced graph of Fig. 3(b). This is where 
k'^ is colhnear to p'^ — see eq. (11). In this region the gluon is moving slowly relative 
to the quark, pj^, and we may regard the gluon and the virtual quark pA + k as being 
given a large boost from the center-of-mass frame. Thus the term in 

{^^ + ^ + M)j\{pa), 

with A = + is by far the largest. 

It follows that the sum over u in (25) is dominated by the term with u = —. 

(Remember that in light-cone coordinates the metric is non-diagonal: = = 

1, 9++ = 9— = 0.) 

Note that we cannot say that the sum over A is dominated by A = +, because 
of the kxki, terms in the gluon propagator. (Even if we set ^ = 1, so that we used 
Feynman gauge, such terms would arise when we consider graphs with vacuum 
polarization corrections for the gluon). 

A more general argument giving the same answer can be made by treating 
all the 7- matrices as order 1. We wish to see how the numerator terms with large 
components, viz. p'\, k~^, p^ contribute. To do this, we anticommute all 7~'s to the 
left and all 7"'"'s to the right and use 7"*"^ = 7"^ = 0. Then we use the mass-shell 
conditions 

i^^-M)uipA) = 0, ^^^^ 
vipB)i^B + M) = 0. 

Since k~^ and p\ always multiply a 7" in the numerator and since p^ always mul- 
tiplies a 7"'", we find that the large terms only arise from anticommuting a 7" with 
the 7'^ or a 7"'" with the 7^^. 

Let us write the numerator as 

Yi-h + ^ + M)Y{h + ^ + ^h'N^u = (28) 

with B'^ = 7'^(— j^R + ^ + M). To simplify this, we use a beautiful trick formalized by 
Grammer and Yennie^. It starts by making the following string of approximations 

B^A^ ~ B-A+ 

= ■^B-k+A+ 

k+ 

1 (29) 

k+ 

= k.Bp^. 
k ■ Ub 
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Here we rely on the facts that + components of A'^ and k'^ are their largest and that 
the — component of B" is not much smaller than its other components. To put the 
result in covariant form, we have defined to be a light-like vector with = 1, 
Ub=uI = 0. 

We now have a factor k^, times the lower vertex. This is a standard situation 
where Ward identities can be used. In the present case the result is easy to derive: 

=v [{-^^ + ]/,^M) + {i>^ + M)] i-^^ + ^ + M) (30) 
=v [{PB - kf - M2] , 

where we used the on-shell condition (27) for the wave function. The factor {pb — 
kY — cancels the antiquark propagator and we find that if k^ is restricted to be 
coUinear to p^, then 



\2 



167r4 Jk coUinear ^ . _ ^2) [^2 _ (^^ + ^) 

with errors being smaller by a power of l/Q (or A, where A is the small scale factor 
in eq. (11)). 

The coupling of the gluon to the antiquark has become featureless; it is in 
fact insensitive to the spin and energy of the antiquark. All the gluon sees is the 
direction and charge of the antiquark. We do not need an ie prescription for the 
pole of 1/k ■ ub at k~^ = 0, for k~^ is always large in the collinear region. 

An exactly similar result holds for the opposite collinear region (Fig. 3(c)) 



~ 167r4 colUnear - {p^ - kY]{m? - k^ - ie){uA ' k) ' ' ^ ' 



Furthermore a slightly different result holds if k^ is in the soft region, sym- 
bolized by the reduced graph Fig. 3(d): 

r~^/ d^k, — ;^M^^^r^>(PA)Ar,.(A:) ^_ 

IQti^ Jsoft [Ub ■ k — ie) [m^ — k^ — ie) [ua ■ k + ie) 

The only subtlety in the derivation of this equation is that we must assume that 
all components of k are comparable (or at least that \k~^/k~\, \k~/k~^\ <^ Q'^/M'^, 
\k~^k~\ > k^). However, there is a leading contribution when k~^ and k~ are of 
order X'^Q and kx is of order XQ, with A a small quantity. This is the Glauber 
regiorPSl^ and there none of the approximations (31) to (33) is valid. Now, in the 
Glauber region — fc^ ~ + k^, independently of k~^ and k~ . So we can get out 
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of this region by deforming the contours of integration over fc"*" and k~ away from 
the poles in the quark and antiquark propagators to where at least one of (31) to 
(33) is valid. To indicate the direction of deformation, we introduced the ze's with 
the ua ■ k and ub ■ k denominators. 

Note that there are overlap regions for k where two or more of the approxi- 
mations (31) to (33) are simultaneously valid. We will see shortly how to avoid the 
double counting that this could give in the factorization theorem. 

The physics of the general case is visible from the one- loop case. The simpli- 
fication in going from eq. (25) to any of (31) to (33) is to replace one or both quark 
lines by an eikonal approximation. That is, the quark is replaced by a source of the 
appropriate charge that exists along a light-like line in either the -|- or — direction, 
and recoil of the approximated quark is neglected. What is happening is that there 
is a large relative rapidity between the gluon and the approximated quark. The 
gluon only sees a Lorentz-contracted object of a certain charge moving at the speed 
of light in a certain direction. On the other hand the quark only sees the gluon for 
only a short time in the quark's rest frame immediately before the annihilation. 

5.3 Leading regions for general graph 

The manipulations in the preceding sections have succeeded in simplifying the in- 
tegrand of (25) in the regions where fc^ remains close to mass-shell as Q goes to 
infinity. We will use these results, and their generalization to higher order to con- 
struct a useful factorized form for the complete form factor. 

The first step (in the general case) is to see that, for a general graph F, all the 
regions that give leading-twist contributions hav^^ the form of Fig. 7. Each region 
is specified by dividing the graph into four subgraphs, designated "jet- A", "jet-B", 
"soft", and "hard". In the subgraphs jet-A and jet-B, the momenta of the internal 
lines satisfy \k^\ 3> \k~\ and \k~\ ^ \k'^\ respectively. A quark line from jet-A and 
an antiquark from jet-B enter the hard subgraph, together with arbitrarily many 
gluons. The hard subgraph also includes the vertex for the current. The momenta 
of the internal lines of the hard subgraph satisfy 3> M^. The soft subgraph 
consists of lines all of whose momentum components are much less than Q. The 
external lines of the soft subgraph are all gluons and attach to one or other of the 
jet subgraphs. Some regions have no soft subgraph. 

Each of the jet subgraphs and the hard subgraph is connected. The soft 
subgraph, if present, may consist of more than one connected component; but each 
of its components must be joined to both jet-A and jet-B. 

Note that within the jet and soft subgraphs there may be loops with large 
ultra-violet momenta. These make up vertices of the reduced graphs (as in Fig. 1). 
These reduced vertices are of the same form as the ordinary vertices of the theory if 
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Figure 7: Leading regions for quark form factor. 



the corresponding regions of momentum space are to give leading twist contributions. 

We must be more precise about the regions represented by Fig. 7. The 
momenta satisfy the following requirements. The value of\k~^ /k^\ for the momentum 
of a soft line must be much less than the value of this ratio for the lines in the jet-A 
subgraph and must be much greater than the ratio for the lines in the jct-i? subgraph. 
This ensures that the Grammcr-Yennie approximation is applicable to the coupling 
of the soft lines to the coUinear lines. The momenta in the hard subgraph must have 
virtualities that are much greater than for the momenta in the soft and coUinear 
subgraphs. This ensures that the Grammer-Yennie approximation applies to the 
coupling of coUinear gluons to the hard part. 

The next step in the proof is to use the same approximations of the Grammer- 
Yennie type that we used for the one-loop graph. After use of Ward identities, this 
will give a factorization. Then we will write the factors in terms of matrix elements 
of certain operators. In general, a given part of the space of loop momenta may be 
in the intersection of several different regions of the form of Fig. 7. We will have 
to make an arbitrary choice of which region to use. The resulting factorization will 
involve graphs with momenta restricted to certain regions of momentum space. We 
will convert this intermediate factorization to a more useful factorization by showing 
that operator formulae representing the second factorization can be converted to the 
same intermediate factorization. 

The resulting factorization will still not be in a form that allows perturbative 
calculations without large logarithms. But it will enable us to derive a differential 
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Figure 8: Grammer-Yennie approximation and elementary Ward identity. 

equation for the Q-dependence of the form factor. The solution of this equation will 
be our ultimate result, and will have all the logarithms separated out. 



5.4 Factorization 



When a gluon collinear to A attaches to the hard subgraph, we use approximation 
(29) where now B'^ denotes the elementary vertex where the collinear gluon attaches 
to the hard subgraph. We then use a Ward identity of the sort illustrated in Fig. 8. 
We let the vertex be B'^ — igY and the momenta of the quarks be I and I + k. Then 



k ■ Ub 



i{l + M) 



k ■ Ub 



[^ + /-M-(/-M)] 



i(/ + M) 



i{l + M) 
k-UB P-M^ 



k ■ Ub 



i{]/^ + l + M) 
{k + lf-M^ 



(34) 

On the right of Fig. 8 we use the double line to denote the eikonal propagator i/k-UB- 
When we sum over all ways of attaching the jet gluons to the hard subgraph, there 
is a whole set of cancellations and the effect is to take the gluons to the A side of 
the hard subgraph, as depicted in Fig. 9. In obtaining this we have used eq. (34) 
repeatedly. Each gluon has a factor i/k ■ ub- Then, we use identities like 



ki ■ Ub 



k2 ■ Ub 



{ki + k2) ■ Ub 



ki ■ Ub 



+ 



{ki + k2) ■ Ub 



k2 ■ Ub 



(35) 



to write the gluon attachments as if they are to a single line with an eikonal prop- 
agator. Physically, Fig. 9 is telling us that the gluons collinear to A only see the 
charge and direction of the antiquark. 

Similar arguments applied to the attachments of the gluons collinear to B 
to the hard part and to the attachments of the soft gluons to the jets give Fig. 10. 
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Figure 9: Result of applying Ward identities to attachments of gluons from jet A to 
the ultraviolet subgraph in Fig. 7. 
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Figure 10: Basic factorization of form factor. 



Diagrammatically, Fig. 10 represents a factorization, but with the momenta in the 
subgraphs restricted to particular regions. Notice that since the momenta in the soft 
part are much less than the momenta in the hard part, we ignore the dependence of 
the hard part on the loop momenta that couple the soft to the hard part. 

Our aim now is to construct a formula that exhibits the factorization of 
Fig. 10, that has exphcit operator definitions of the factors, and that has no restric- 
tions on the momenta of the lines in the Feynman graphs. First let us observe that, 
for example, the Feynman rules for both the jet-A subgraph and the cikonal line 
attached to it can be derived from the following matrix element of the quark field 
with a path-ordered exponential of the gluon field: 







Texp 



^0 



dzuB ■ A{—ubz) 



q{0) 



Pa 



(36) 



Now consider the following quantity: 



Ja 



Pa ■ 



;m,M,g,fx] = 



(0|Texp [igfi^J^dz n ■ A{nz)] q{0)\pA) 
(OlTexp [ig/i^ dz n ■ A{nz)] |0) 

X UV renormalization factor. 



(37) 
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Figure 11: Typical graphs with eikonal self interactions. 




(a) 



fbl 



Figure 12: One loop graphs for Ja- 



The numerator is the same as (36) , except that we have replaced by a space- like 
vector = u'^ — u'^. This has the effect of suppressing the contribution of momenta 
for which \k"\ ^ l^i^^l, i-e., momenta collinear to B. The need for the vector n'' 
to be a space-like rather than timc-likc will appear later. The denominator in 
(37) is necessary to cancel graphs, like Fig. 11, with eikonal self-interactions; these 
do not appear in Fig. 10. Finally, since we have removed all restrictions on loop 
momenta, there are ultra-violet divergences in graphs like Fig. 12; these we define 
to be cancelled by renormalization counterterms. 

The one-loop contributions to J a are given in Fig. 12, so that 



Ja = zI^'^u{pa 



|4-2e 



IGtt^' 

-|- UV counterterm -|- 0{al). 



k 



Y{i>A + ^ + M)N^,{k)n-u{pA) 
{m? - A;2) [M2 - (p^ + kf] {n ■ k + ie) 



(38) 

This reproduces the contribution to Fig. 6 of the region where k^ is collinear to A. 
In (38), Z2 is the residue of the pole of the quark propagator. 
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A jet-B factor may be defined similarly: 

7 ( P^ • M ) = (0|rg(0)exp[iff/x^/o°°dz n- A(-nz)] \Pb) 

Jny , ,m,g,i,j- {0\T eMig^^^ J.^ dz n ■ A{-nz)]\0) (39) 

X UV renormalization factor. 

Wc now apply the same argument to Ja and Jb as the one we applied to 
obtain Fig. 10 from Fig. 7. The result has the form 

Ja — (Jet-74) x soft factor x hard factor, 

(40) 

Jb = (Jet-B) x soft factor x hard factor, 

where "Jet-A" and "Jet-S" are the same quantities as in Fig. 10, but the soft and 
hard factors are different. Hence the form factor can be written 

F{Q) = Ja (^^) X Jb (^^) x soft x hard + 0{1/Q^). (41) 



We next recognize that the soft factor in Fig. 10 has the Feynman rules for 







Texp 



/ dz ua- A (—zu'^) exp —igji'^ / dz ub ■ A {—zub) 
Jo J L io 



0), (42) 



with the momenta restricted to the soft region. Suppose we were to define a soft 
factor as exactly (42) without any restrictions on loop momenta. Then there would 
be divergences from regions where gluons become collinear to ua or ub- These 
divergences are caused by the fact that -u^ and Ub are light-like and effectively 
represent an incoming quark and antiquark of infinitely high energy. An example is 
given by the one loop graph, Fig. 13: 

''^ ^d^fe ,,_, ...V!L. (43) 



(27r)4 J [k- + ie) (m? - k^) {k+ - ie) 

The collinear divergences come from the regions where k^ /k~ — or cxd with k^k" 
fixed. They are evidently artificial divergences. The actual collinear regions of the 
original form factor have already been taken into account by the factors J a and Jb in 
eq. (41), so we would also be guilty of double-counting if we were to keep exactly (42) 
as our definition of the soft factor. (There is also an ordinary ultraviolet divergence 
from the region where \k'-^\ — > 00; we will deal with this separately.) 

What we must do is to observe that the argument that led to the factorization 
of Fig. 10 for the form factor can also be applied to (42). The result is that the 
collinear parts factorize. This is shown in Fig. 14, where the soft factor is identical 
to the soft factor in eq. (41). Therefore we can write the original form factor as 

F = Ja{pa ■ r?ln^)JB{VB ■ n'/n^)S{M, m, g, //) x hard + 0{1/Q^), (44) 
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Figure 13: One loop graph for eq. 



(42). 




Figure 15: Definition of soft factor. 



where S is the quantity (42) divided by its coUinear divergences: 



S= (0 



Texp 





[(» 


Texp 






Texp 



poo poo 

gjf I dz ua- Ai—zu'X) exp —ig/J," / dz ub ■ A(—zuq) 
Jo J L ^0 

roo "I r roD 1 \ ^ —1 

ig/j," J dz ua- A. {—zua) exp igjj,'' J dz n ■ A (zn) 0)> 

poo "I r POO "I \ ^ 

J dzn-A[—zn) exp —igjfj dzuB-A[—zuB) Oy | 



X (0 



Texp 



igU 
ig^i 



dz n • A izn] 



exp 

o)/o 



-1 



Texp 



igii 



dz n ■ A i—zn) 



. 



(45) 

This is represented by Fig. 15, and the last fine in the equation has as its sole purpose 
the canceling of eikonal self-energies. In (45) we implicitly define UV divergences to 
be cancelled by renormalization counterterms. 

We finally observe that when the form factor F is divided by JaJbS, all its 
coUinear and soft regions have been cancelled. Hence we can define the hard factor 

by 

H{Q,g,fi) = ^mF/{JAJBS). (46) 

m-»0 

Therefore we obtain the factorization 
F{Q,m,M,g,i^) = 

Ja(-QV4; M, m, g, fx) Jb{-QV^] M, m, g, //) S{m, M, g, /x) H{Q; g, /x) (47) 

+ o(i/g2). 

Here Ja, Jb, S and H are defined by Eqs. (37), (39), (45) and (46). 
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We cannot directly use this equation to control the Q- dependence of F. 
However, the dependence of J a and Jb on Q is through the vector n'^, since 
= —ApA ■ n'^/n^ = —^Pb ■ v? jr?. What we will do in the next section is to 
compute the Q-dependence of J a and Jb by differentiating each with respect to n'^, 
holding the physical momenta p\ and fixed. The Q-dependence of H will be 
renormalization-group controlled just as for the Q-dependent factor in 0^ theory. 

The reason for using a spacelike vector can now be explained. On the 
various occasions that we factor out a soft region we needed to deform integrals 
over gluon momenta away from the Glauber region \k'^'k~\ <C just as in deriving 
eq. (33). For consistency all these deformations must be in the same direction. In 
the case of a soft region for the jet factors J a and Jb this means that the n ■ k 
denominator must be {k~ — k^ + ie) to give the same direction of deformation as in 



6 EVOLUTION EQUATION 

In this section we will derive equations for the Q-dependence of the form factors, 
first in (0^)6 theory and then in a gauge theory. 



The Q-dependence of the form factor (17) for cf^ theory is under renormalization- 
group control — see eq. (21) or (22). That is, after the factorization is obtained, 
we may change /i to different values in the two factors to eliminate all the large 
logarithms in their perturbation expansions. So we need not bother to derive an 
explicit equation for its Q-dependence. However, for the sake of the comparison 
with the case of a gauge theory, we will nevertheless do so. 

From eq. (22) we have 



The right-hand side of this equation can be expanded in powers of g. The expansion 
contains no logarithms of Q or of the masses, so it is valid to approximate it by the 
first term or two, provided only that g{Q) is small. Since 0^ theory is asymptotically 
free , su ch an approximation is valid for all large enough values of Q. Indeed, given 



(33). 



6.1 



dlnF 
d\nQ 



7 igic2Q)) + P {g{c2Q)) ^ In F (l/c2, gic^Q)) + 0{l/Q^) 
7(^7(c2Q),C2) + 0(l/g2). 



(48) 



thafii 




(49) 
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we find that 

dlnF _ -4/3 

~ In (gvA2) ^ ■ ■ ■ ■ 

Here we have used the value of 7 that can be extracted from the one-loop form 
factor (16): 

7 = Q^eq. (16) + 0(^^) 

9Q (51) 

= -//(647r3)+0(^^). 

Note that the calculation of 7 from the vertex graphs may be done entirely in the 
massless theory. This results in a saving of calculational efi^ort. 

Since the right-hand- side of (48) is of order g'^{Q) with no extra logarithms 
of Q, the variation of the form factor for, say, a doubling of Q is small, of order 
1/ ln((5/A). This is comparable to the scaling violations in deep-inelastic scattering. 
If we integrate over a wide range of Q, say from a value Qq to a value of order Ql/A, 
effects of order unity arise. If both Q/A and Qq/A. are large, then 



F{Q) = F(Qo) 



In(gVA^) 



-2/3 



[1 + O(lnln(go)/In(go))]- (52) 



These results will be used as a standard of comparison when we have derived the 
corresponding results for a gauge theory. 

The accuracy of the results may be systematically improved by calculating 
higher orders in perturbation theory. 



6.2 Gauge theory jet factors 

In a gauge theory, the form factor satisfies the factorization equation (47) . The hard 
(or ultra-violet) factor H has Q-dependence that is renormalization group controlled, 
just as in (0^)6 theory. But there is further Q-dependence in the jet factors. This 
comes essentially from the possibility of emission of gluons of moderate transverse 
momentum in a range of rapidity bounded by the two incoming particles. These 
gluons are divided into what we may term left-movers and right-movers by the vector 
n^. The remaining contributions are put into the soft factor S and the hard factor 
H. The details of the cut-off on rapidity given by the vector n'^ are incorrect in the 
central region of finite center-of-mass rapidity. The errors are compensated by the 
Q-independent factor S if they correspond to quanta of low transverse momentum. 
Quanta of large transverse momentum are also included in all of the factors J a, Jb 
and S, but they are incorrectly approximated; the hard factor H was defined to 
cancel these errors. 

If we were to directly investigate the Q-dependence of the form factor, we 
would have to find its dependence on and p^. We would have to trace the fiow 
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of these extern al mom enta inside Feynman graphs, and this would be a hard task. 
(See Sen's worl^^^^ for details on how to do this.) But if, instead, we examine 
the factors J a and Jb, we see that their Q-dependence comes from dependence on 
Pa ■ T? It? and "Pb ■ /n^ respectively. So it is sufficient to find their dependence 
on n^. This is very much easier because it involves the result of differentiating in 
their definitions the path-ordered exponentials of the gluon field, with respect to 
direction. This gives a very simple result, as we will now see. 

We have for J a 

OJa . ,ydJA , . 

where 6n^ is a backward-pointing time-like vector, normalized so that = —rP 
and n ■ 6n = 0. With our previous representation, where n'^ = u'^ — u'^, we have 
6n'^ = — — M^. The light-like vectors and are defined hj u\ = u]^ = 1, 
UA = ut = u\ = u^B = 0- 

The Q-dependence of Jb is obtained by the opposite variation of n^: 

Since this will result in an equation identical to that for the Q-dependence of Ja, 
we will restrict our attention to Ja- 

The result of differentiating the path-ordered exponential with respect to its 
direction is 



oo 



igif I dz n ■ A{nz) 

' (55) 

/■oo r /-oo n 

/ dz [5n ■ A + z5n ■ dA ■ n] x exp ign"^ / dz n ■ A , 
Jo L "'0 J 



where we used the canonical equal-time commutation relations of the gluon field 
to commute A^ with A'^ and dA^/dt with A^. The resulting Feynman rules for 
dJA/d\n.Q are exhibited in Fig. 16 and Fig. 17. The version shown in the second 
line of Fig. 16 comes from the second derivation which we now give. 

Suppose in a given Feynman graph (Fig. 18(a)) for J a there are gluons 
attaching to the eikonal line. We have a factor 

—gn^^ —gn^'^ —gn^^'^ 

(56) 



{ki ■ n + ie) {ki ■ n + k2 ■ n + ie) {ki ■ n + k2 ■ n + . . . k^ ■ n + ie) 

Now let us sum over all graphs which are the same as the first one except for having 
the gluons permuted. The result is to replace (56) by 



-gn^^ —gn^'^ —gn^^ 



{ki ■ n + ie) {k2 ■ n + ie) {k^ ■ n + ie) 
26 



(57) 
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ign^^ -, ^- ]■■■[ ign'"^ ■ 



ki ■ n + iej V k\ ■ n-\- . . .kN ■ n -\- ie 
Figure 17: Feynman rules for Fig. 16. 



Figure 18: Sum over attachments to eikonal line gives eq. (57). 

and is illustrated in Fig. 18(b). Finally, we take the derivative given on the right 
of eq. (53). This results in a derivative for each of the N gluons, as symbolized in 
Fig. 16, each of the derivatives having the form 

Jll^ = Sn''^(^P^\^- ^—-{5n-kn^-n-k5n^), (58) 

n-k + ie dn^Kn-k + ieJ (k-n + ie)^^ ' ^ ' 

as stated in the Feynman rules. Fig. 17. 

To derive an evolution equation we now use the same style of argument that 
we used to derive our first factorization (47). We can again divide momenta into 
collinear, soft and hard. (The collinear momenta here are only those collinear to A.) 
The simplification that now occurs is that a collinear momentum cannot enter the 
differentiated vertex. This is easy to see, since from (58) we find a factor 

which is suppressed by a large factor compared with 

A-n^A--A^ (60) 
from a regular eikonal vertex. Here k^ and are vectors collinear to p^. 
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Figure 19: Leading regions for Ja- 




Figure 20: Factor obtained from Fig. 19 when k is (a) soft, or (b) hard. 

Therefore, to leading twist, only soft and hard momenta attach to the dif- 
ferentiated vertex. The result is Fig. 19, which is analogous to Fig. 7 for the form 
factor. 

Whenever /c^ is soft we can use the Grammer-Yennie method to obtain the 
factor shown in Fig. 20(a). Whenever is hard its line disappears into a hard 
subgraph; then we apply the derivation of Fig. 9 to Fig. 20(b). This gives 

dJ 

— — — = Ja X [soft + hard] + higher twist. (61) 



6.3 Operator form of jet evolution 

It will be convenient to rewrite the soft factor in (61) so that it has an explicit 
definition in terms of a matrix element of an operator, just as we rewrote the jet 
and soft factors in Fig. 10 to give eq. (47). The same method of argument as for 
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Figure 21: Definition of K. 



Fig. 10 gives 



dlnJA 1 



K (m, M, ji, g{iJi)) + hard. 



(62) 



d\uQ 2" 

We have defined the soft term K to have a factor 1/2, since it will be multiplied by 
2 when we write the evolution equation for the form factor F. Furthermore, this 
equation will have an extra hard term coming from the factor H in eq. (47), so we 
do not bother to name the hard term in (62). 

The definition of the quantity K is given in Fig. 21. It can be written as 



X = 2(0 



Texp 



igii 



/ dz ua • A (—zua) exp ign" / dz n ■ A(zn) 
Jo J L Jo 



igji" / dz\5n- A-\- z5n- din ■ A)] 

Jo 



T exp 



igU 



dz Ua • A{—zua) 



exp 



igfx 



dz n • A{zn) 



with renormalization. 



(63) 

No restriction is placed on the internal momenta of so we must add UV coun- 
terterms. Sub divergences in K all correspond to the usual renormalizations of the 
interaction, and the remaining divergence is an overall divergence, so we must define 
it to be cancelled by an additive counterterm: 



K{m, M, g, fi) = K^^are + counterterm 
= Ky,^re + SK{g,e). 



(64) 



All soft contributions to din Ja/OIiiQ are contained in K together with some hard 
contributions. 
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6.4 Gauge theory form factor 



We can now write an equation for the form factor 

d\nF dlnJA dlnJe dlnH , , 

dh[Q = + 9Mg + 9h^ + '"^"^ 

In H 

= Kim, M, fif, u) + 2 X hard + -— — — + higher twist 
^ ^' d\aQ 

= K{m, M, g, jj,) + G{Q/fi; g) + higher twist. (65) 

Since the dependence on masses has been separated from the dependence on Q on 
the right-hand side of this equation, we wiU have an effective calculation of the large- 
Q behavior of the form factor once we know the renormalization-group equation for 
K. 

This equation is easy to derive, since i^bare is renormalization group invariant: 
{^dK\yaTe/dfi = 0). So from eq. (64) we find that 

/i— = fi—6K 
d/i d/i 

= -1k{9)- (66) 

The anomalous dimension 7^ is finite at e = 0, and if we use minimal subtraction, 
't Hooft's methodJI^ show that it is independent of e. The anomalous dimension of 
the form factor is zero, since the renormalized operator [qj^q] has zero anomalous 
dimension: 

Thus eq. (65) gives the anomalous dimension of G: 

fi^G{Q/^;g{^)) = - \im 

an m,M^O ajJ, 

= IK (c/(/i)) • 

We can now write the evolution equation in a form with no large logarithms: 
d\nF 



d\nQ 



Q da' 

-TlK {gifi')) - G (I/C2; g{c2Q)) - K (m, M, g{^^),^^) 

M /i 



(69) 



Here we displayed an overall factor —1, because the dominant term is the integral 
over the positive one-loop value of 77^. 
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7 INTERPRETATION 



The physical interpretation of eq. (69) and its derivation is as follows: 

When the total energy Q is increased, the phase space available increases 
for virtual quanta inside the form factor. Let us consider the size of the phase 
space split up into different ranges of transverse momentum. At low transverse 
momentum the range of phase-space is governed by the rapidity between the quark 
and antiquark — this gives the K term. At large transverse momentum, of order 
Q, the quanta can only have finite rapidity, but the range of transverse momentum 
increases with Q — this gives the G-term. Finally, the intermediate range is filled 
in by the anomalous- dimension term. 

The precise form — e.g., the Q-independence of K — and the detailed deriva- 
tion rely on the fact that a particle (virtual or real) can only probe details of another 
particle (e.g. the initial quark or antiquark) if the relative rapidity is low. At large 
relative rapidity there is not sufficient proper time to get a detailed picture. Indeed 
the only elementary particle that can even interact at all across a large rapidity gap 
is the spin-1 gluon. Then it just measures the total charge and the direction of the 
probed particle. A coherent sum over the detailed structure is needed to give this 
result in perturbation theory, the result being formalized in the Ward identities. It 
is the need to sum over a set of Feynman graphs to get the physical answer that 
results in the technical complication of our derivation. 

The key to understanding the derivation is to ask what happens to quanta 
inside the form factor when Q is increased by boosting the quark and antiquark in 
opposite directions. We first examine quanta with transverse momenta much less 
than Q. The argument in the previous paragraph indicates that the interactions of 
these quanta with other quanta of very different rapidity, do not depend on the size 
of the rapidity gap. Thus the part of the change in F that comes from quanta of 
low transverse momenta can be found by measuring the quanta that come into fill 
the interior of the increased rapidity range. 

To make this measurement we first measure the contents of the incoming 
quark down to some finite rapidity, using the operator 



Then we differentiate with respect to n'^, to find the change caused by increasing 
the range of rapidity. 

The remaining part of the variation of the form factor with Q comes from the 
short-distance regime of large transverse momentum. This region is well understood. 
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8 SOLUTIONS AND CALCULATIONS 



It is easy to solve eq. (69), with the result 



F{Q)=F{m,M,g,fi)exp\ - ^ In 




+ K{m,M,g{fi),fi) ln(c2Q//i) 




by the form factor at low Q. If is chosen to be of order the masses then this 
is a form in which no large logarithms appear in the coefficients of perturbation 
expansions; the logarithms in the unimproved perturbation series for F are either 
explicitly in the exponent in (70) or are generated by the integration over /i'. 

The leading logarithmic approximation to F{Q) is obtained by using the 
result (see later) that •jk = fi'^/(27r^) + 0{g^) and by writing g{fi') in terms of the 
running coupling, g{fi), at a fixed scale. Then the highest power of a logarithm of 
Q is obtained from the lowest order term in the exponent. The result is 



In QCD, formulae like (70) can be derived for a number of important cases, 
such as the transverse-momentum distribution of the Drell-Yan procesJ^^. Since 
QCD is asymptotically free, we can do an effective calculation from low orders of 
perturbation theory if Q is large. Non-leading logarithmic corrections are thereby 
tamed. 

8.1 Calculations 

We now consider how to calculate the quantities K, G and, particularly, '-fx that 
appear in the exponent in (70) and on the right-hand side of the evolution equation 
(69). One method is to start from the Feynman graphs for the form factor. Then 
the results of these calculations are compared with the general form eq. (70). This 
determines K and G (and hence jk) except for an ambiguity of adding some function 
of g to K and subtracting the same function from G. This ambiguity is the same as 
the renormalization-scheme ambiguity for K and G, and as such does not directly 
pertain to physics: the physics resides in the dependence of the functions K and 
G on their other arguments. Especially at higher-order this procedure is rather 
lengthy. 



F = exp 




In^ + non-leading logs 



(71) 
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Now the terms in the exponent with the largest powers of logarithms are 
the most important — in particular the 7x term. So a short-cut can be made by 
calculating K. At two- loop order there is only one easy graph for i^T — a vacuum 
polarization correction — but for the form factor there are five harder graphs. The 
anomalous dimension 7^ is obtained from K and in particular from its ultra-violet 
divergence. 

In addition, the form (70) implies many relations between the logarithms of 

Q in different orders. These are not manifest in a direct calculation. However their 
validity provides nontrivial tests of calculations. Nevertheless all the coefficients can 
be obtained by a direct calculation. For example, from an evaluation of the one-loop 
vertex, eq. (25), one can show that 



In — ^ + 0{al) + higher twist. (72) 



This is evidently of the form of the evolution equation, eq. (65). It implies that 

(73) 



TT 



In(m2) + C +0{al), 



TT 



ln(-g^) + C +0{at), 



where C is a constant. The value of C is not a priori fixed, and a change of the 
constant corresponds to a change of renormalization scheme for K. We choose to 
resolve the ambiguity by using MS renormalization applied to a direct calculation 
of K from its Feynman rules. Given that to 1-loop 

Stt^ ^ J {m? -k^ -ie){uA-k + ie){n-k + ieY 

= J ^""^ im^-k^-^e)i-k^ + ^er + 

_ (74) 

it is fairly easy to perform the integrals. The result, with MS renormalization, is 
that 

K^^]n(myf/)+0{a^,), (75) 

from which follows 

7K = 2^ + 0(a,^). (76) 

TT 

It is left as an exercise for the reader to show that the sole 2-loop graph gives 
the 0{a^) term in jk- 
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There is a lot of information in these results, even without the two-loop result 
(77). For example, let us expand InF in powers of t = ln(— Q^Z/i^) 



\nF = — (Cut^ + Cut + Cio) 



+ 0(1/Q^), 

where the coefficients may depend on m, M and /i, but not on Q. The leading 
logarithm results imply that C24 = 0. Our formula (65) implies considerably more. 
Now, from (78) we have 

"•(4C„« + 2C„) 



d\nQ vr 



+ ( ' f 6C23t' + AC22t + 2C21) ^^^^ 



+ . . . . 

In order that G in eq. (65) be independent of the masses m and M, C12, C23 and 
C22 must be independent of m and M (and hence of /i). Furthermore, once one puts 
in the one-loop values, the requirement that G satisfies its renormalization group 
equation implies that 

C23 = (80) 

Hence the new information for the form factor F at 2 loops is 2 logarithms 
down from the leading logarithm, i.e. it is in C22 and the less leading coefficients, 
C21 and 6*20- The double logarithm coefficient C22 is related to the two- loop term 
in '-yxi given in eq. (77); this was the result of a relatively easy calculation. Hence 

= 1. (81) 

The remaining information, for which a full two-loop calculation of the form factor 
is needed, is in the terms with one and no logarithms of Q. These are three and 
four logarithms down from the leading In"^ Q term. 



8.2 Comparison with other vi^ork 

One can verify (77) and (80) from the calculation of Barbieri et a^^. Note that in 
this calculation one must change renormalization prescription first. 

Korthals-Altes and de RafaeP^ made a conjecture about an evolution equa- 
tion for the form factor we are discussing. Their conjecture is that {Qd/dQ — 
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(3 d/dus) InF is linear in InQ. It can be checked that their conjecture is imphed by 
our eq. (69). However the converse is not true: Their conjectured result does not 
imply our eq. (65) with its specific dependences on masses. 

A number of calculations of comparable quantities in QCD have been made. 
As will be discussed below, generalizations of our formulae apply not only to a simple 
quark form factor but notably also to the transverse momentum distribution in 
the Drell-Yan process and to two-hadron-inclusive production in e^e~ annihilation. 
(This last includes the energy-energy correlation as a special case.) The anomalous 
dimension 7^ is common to all these processes. 

The electromagnetic form factor of a quark in massless QCD also satisfies 
our equation (65) or (69), as shown by SerP^. The coefficients are now 



(82) 

Here n/ is the number of quark flavors, while Tp, Cp and Ca are the usual group 
theory coefficients. (Tp = |, Cf = | and Ca = 3 for QCD, while Tp = I, Cp = I 
and Ca = for QED.) Since there are no masses, -ft^bare is zero, and the renormalized 
K equals its MS counterterm. The resulting pole, as displayed in eq. (82), represents 
the infra-red divergence in K. Those terms that appear in the abelian case are given 
by our earlier calculations. The only purely nonabelian term in the order to which 
we work in eq. (82) is the two- loop CpCa term. We have deduced its value from 
I 3U I 31 I 32 |^ J explain later. 

The calculation in QCD that can most directly be compared with eq. (82) 
is by GonsalveJ^ who has calculated precisely the quark form factor in QCD, at 
two-loop order. (The purpose of doing this is that one can use the deduced value 
of in other processes.) Gonsalves' results do not obey the correct evolution 
equation, which should hold in QCD as well in QED. Note that his renormalization 
prescription differs in detail from both MS and MS. The agreement between the 
other calculations indicates that there must be an error in Gonsalves' calculation. 



K 



— Ljp 

TT 



7 -|- ln(47r 

e 



C = -Cp 

71 

7x = 2Cp^ + 

TT 



In {Qy^^^) - I + O («: 



8.3 Infrared divergences in QCD 



Korchemskii and Radyushkin ''^' I have studied the infrared divergences of the elec- 
tromagnetic form factor of a quark in QCD at large Q. (Their ultimate ainJ^ is to 
study the full Sudakov problem in QCD.) Consequently their methods have much in 
common with the work described in this paper. Indeed their results are written in 
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terms of path-ordered exponentials that are similar to the ones used in this article. 
Moreover they have derived the necessary generalization of the Grammer-Yennie 
method to the nonabelian case. 

In an abelian theory, the infrared divergences (as the gluon mass goes to 
zero) are rather simple: they form a factor which is the exponential of the one loop 
infrared divergence. In eq. (69), for the Q-dependence of the form factor, all the 
infrared divergences are in the term K. Gluon self couplings in an abelian theory 

are induced solely by quark loops, and these loops are suppressed when the gluon 
momenta go to zero and the quark mass is nonzero. Hence only the one loop part 
of 7ii- is needed for the calculation of the infrared divergences in the abelian theory. 

In QCD the infrared divergences are much more complicated, since the 
gluon self couplings are not so suppressed at zero momentum. Korchemskii and 
Radyushkin show that the infrared divergences form a factor: 

_ roc roo 

Mso{t = {0\TPexp[-igJ^ ds p'XA^ipAs)] P exp[ig ds p%A^ipBs)]\0)m. (83) 

The subscript TR' means that integrations are restricted to the infra-red region. 
When = {pA —PbY gets large, the infra-red behavior is governed by the anoma- 
lous dimension of the cusp, Fcusp- The ability to do systematic perturbative calcu- 
lations at large Q relies on the property, proved by Korchemskii and Radyushkin, 
that Fcusp is linear in ln((5) for large Q: 

Teusp = A{a,) HQ/M) + B{a,) + 0{l/Q^). (84) 

This linearity is implied by our results (if it is assumed that they extend to QCD). 
Indeed A{as) in eq. (84) is the same as our 7^^:.^ The proof is simple: 

d _ d d , 

2 ^"«P - pi^^n2A^^.. -^^"f* 



aing2 ^"^P ainQMln// 

d d 
dlnjidlnQ'^ 



\yiM 



soft 



dln/x 



dln^u 



K 



Here we have employed the factorization theorem for In Alsoft that is analogous to 
the one for the quark form factor. The soft term K is the same as for the form 



^Note that in going to the regime of infrared divergences, it is necessary to compute the anoma- 
lous dimensions in the effective low-energy theory that exhibits the decoupling of massive quarks. 
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factor. But the hard term Gkr is different. Indeed, since it is a pure ultraviolet 
quantity, with no scale dependence, it is zero. 



This result enables us to obtain the non-abelian part of at two loop order; 
this is the term in eq. (82) that is proportional to CfCa- Since Korchemskii and 
Radyushkin investigate the infrared divergences of the form factor with massive 
quarks, they cannot calculate the CpTp term. 

8.4 Other QCD calculations 

Kodaira and TrentaduJ^ considered the energy-energy correlation in e"'"e~ annihi- 
lation. They worked with a different formalism for the Sudakov form factor. Their 
calculation was the first from which the nonabelian part of the two-loop value for 
can be deduced. They also agree with the abelian part, as calculated directly 
from K. 

Davies and Stirlin^^ have calculated the Drell-Yan cross section at order 
Og. They deduce 7i<- and the equivalent of K and G. They confirm the value for -^k 
given in eq. (82). 

9 APPLICATIONS TO QCD 

As has already been noted, there are many cases in QCD where something like a 
Sudakov form factor enters. The most straightforward extension of the results in 
this article is to transverse momentum distributions. The transverse-momentum qx 
is a third important scale for the cross-section, in addition to the total energy Q and 
the hadron mass-scale. Two logarithms of Q/qx per loop are present in Feynman 
graphs. 

In the case of two-particle inclusive cross-sections in e^e~ annihilation, Collins 
and SopeJ^ derived an equation generalizing Eqs. (65) and (69). The same anoma- 
lous dimension '-y^ niakes its appearance. Technically the main difference between 
^ and the treatment in the present article, aside from having a non-abelian gauge 
group, was that there we used an axial gauge n ■ A = instead of Feynman gauge. 
This resulted in a nice simplification. For example, in eq. (37), the line integral of 
the gluon field is zero, so that in^ the definition of J a would have been 



with the dependence on n'^ now being a dependence on the choice of gauge. 

My treatment of the Sudakov form factor in^^ used Coulomb gauge, which 
behaves for this purpose rather like the axial gauge. In either gauge, explicit Feyn- 
man graph calculations are made more difficult than in covariant gauge by the 




(86) 
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complicated form of the numerator of the gluon propagator. In axial gauge we have 



where 'PV denotes the principal- value prescription for the singularities at n - k = 0. 
In Coulomb gauge, we have 



where 6n'^ is the vector defined just below eq. (53). 

However, there are more fundamental disadvantages than calculational com- 
plexity to use of these physical gauges. It is very hard to define higher order graphs 
in the axial gauge because of the need to multiply principal values. To overcome 
this, considerable complication in the Feynman rules is necessarjJ^, and the sim- 
plicity of the Ward identities is no longer clear. There are also complications in the 
Feynman rules in Coulomb gauge beyond 2-loop ordeJ^. In both cases, it is not 
clear that a complete and correct all-orders derivation can be given easily. 

We expect corresponding results to the ones for the energy-energy correlation 
to hold for the Drell-Yan process; they have been formulated by Collins, Soper and 
StermarP. This work, because it entails a complete treatment of a factorization 
theorem for a process at low transverse momentum, includes treatment of intrinsic 
transverse momentum within QCD. In other work of that period, based on leading 
logarithmic formulations, intrinsic transverse momentum tends to appear as an ad 
hoc phenomenological modification to the basic formula for the cross section. 

Davies and StirlinJ^^^ have applied this formalism phenomenologically. 
Altarelli et a/.^ have also performed phenomenological calculations, but without 
the full treatment of the intrinsic transverse momentum effects. 

A further disadvantage to using the physical gauges appear^ when we try 
to derive results for the Drell-Yan cross-section. The problem is that singularities 
in the numerators of the gluon propagators (87) or (88) wreck the derivation of the 
form of the leading regions. Specifically, we need contour-deformation arguments 
generalizing those which we summarized at the end of Sec. 5, and these are invalid in 
an axial or Coulomb gauge. So Collins, Soper and Sterman^ were forced to the use 
of a covariant gauge, at the price of some extra technicalities in the proofs. At the 
same time, the proofs come out to be cleaner. It is a generalization of the method 
of ^ that is used in the present article. 

Another line of development comes from realizing that similar physical phe- 
nomena to those in the Sudakov form factor occur inside of amplitudes for scattering 
in the Regge region. SeiP has produced very important results in this He 
used Coulomb gauge. 
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